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ABSTRACT 



A high-resolution, multi-level, primitive equation ocean model is used to examine 
the response of an idealized, flat-bottomed, eastern boundary oceanic regime on a beta- 
plane to both steady and daily-varying climatological wind forcing. The area of study 
is a coastal region within the Chile Current System from 22® S to 34® S. When steady 
wind forcing is used, an equatorw'ard surface current and poleward flowing undercurrent 
develop. Eddies are also generated, with initial formation in the poleward end of the 
domain. When daily-varying wind forcing is used, there is large spatial variability in the 
oceanic response. A relatively w^eak polew^ard flowing undercurrent appears, first in the 
polew'ard part of the domain. An equatorward surface current also develops and inten- 
sifies during the upwelling season. Eddies are generated and develop farther polew'ard 
in the domain than in the first experiment. The eddies are largest in the equatorward 
end of the domain. The eddy motion is closely tied to seasonal influences, with cyclonic 
(divergent) eddies traversing towards areas of higher dynamic heights when there are 
divergent wind fields present, and vice versa. In addition to the currents and eddies, 
upwelling and cold filaments are evident in both experiments. 
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I. INTRODUCTION 



The eastern boundary current system adjacent to the west coast of Chile, although 
well known for its biological productivity, has (until recently) attracted relatively little 
investigation into the physical oceanography of the region, despite an obvious relation- 
ship between the fertile waters in the region, and the fluid dynamics that make them so. 
The mean climatological current system is comprised of several major features, includ- 
ing: the equatorward flowing surface current called the Chile-Peru Current, and a 

poleward flowing countercurrent and/or undercurrent known as both the Gunther Cur- 
rent and the Peru-Chile Undercurrent. 

Eddies, jets and other instabilities that form along this and other eastern boundary 
current systems (e.g., the California Current System) have begun to generate a signif- 
icant amount of interest. Process-oriented studies using numerical models, such as those 
spearheaded by Batteen et al. (1989), attempt to quantify the forces that create these 
features. 

Because the Chile Current System has not been previously investigated with nu- 
merical models, the objective of this thesis will be to investigate the effects that wind 
forcing has on the generation of currents, eddies and filaments in the Chile Current sys- 
tem. A high-resolution, multi-level, primitive equation model will be used to simulate 
the effects of wind stress on the eastern boundary' current region off Chile. The numer- 
ical model used for this study is one developed by Haney (1985), and adapted for eastern 
boundary current regions by Batteen (1989). 

The study is organized as follows. In Chapter II we review’ the oceanic and wind 
regime off the coast of Chile. In Chapter III we discuss the numerical model used, as 
well as the specific initial and experimental conditions. In Chapter IV we describe the 
analysis techniques used to investigate the dynamical reasons for eddy generation. The 
modeling results are presented in Chapter V, and a summary is given in Chapter VI. 
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II. BACKGROUND 



A. REGIONAL DESCRIPTION 
1 . Surface current 

The mean climatological ocean adjacent to the coast of Chile is an eastern 
boundary' current regime that consists of a slow-moving (~ 10 cmjs) surface current 
overlying a relatively strong (~ 25 cm/5) and steady coastal undercurrent (Wooster and 
Reid, 1963). The surface current is similar to other eastern boundary currents in that it 
is a sluggish, wide and generally shallow stream moving equatorw’ard along a west coast 
(Wooster and Reid, 1963). Mass transports have been estimated to be between 
17.4 X 10* - 19.0 X 10* m/5 (Wooster and Reid, 1963). The primary source of the Chile- 
Peru Current is, as expected, the West Wind Drift, because the South American conti- 
nent extends far enough south to intercept at least part of the flow, and divert it towards 
the north. Because the high latitudes are areas where precipitation typically exceeds 
evaporation, the northward flowing waters are relatively low in salinity, averaging less 
than 34.5 ppt. As the flow moves towards the north, salinity levels increase because of 
mixing and the general increase in evaporation along the way (Wooster and Reid, 1963). 

Tchernia (1980) and Gunther (1936) characterize the Chile Current System as 
consisting of two equatorw'ard flowing currents, separated by a seasonal surface 
countercurrent. The outer current is usually referred to by those who adhere to this "two 
current" philosophy as the Oceanic Chile-Peru Current; it extends to depths of never 
more than 1000 m (Wooster and Reid, 1963), and usually no more than 700 m (Tchernia 
1980). The inner surface current is known as the Coastal Chile-Peru Current, with a 
maximum depth of ~ 200 m (Tchernia, 1980). 

Another description of the Chile-Peru Current is that it is a single, broad, slow, 
and sometimes intermittent surface current overlying a poleward flowing undercurrent 
(see Fonseca, 1989). In addition, Wooster and Reid (1963) describe the Chile-Peru 
Current as the principal component of the meridional circulation in the South Pacific. 

Yet another view of the Chile Current System is a combination of the previous 
di.scussions. Fonseca (1989) describes the current as being primarily one body of water; 
however, he states that in the southern domain the current is physically split into two 
components, due to either a countercurrent as described by Tchernia (1980), or to an 
undercurrent that has broached the surface (discussed in the next section). 
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2. Peru-Chile Undercurrent 

Beneath the surface flow is a steady, poleward undercurrent, generally known 
as the Peru-Chile Undercurrent (PCU). The undercurrent is strong and steady enough 
to be considered by some (e.g., Fonseca, 1989) to be the dominant ocean circulation 
feature in the Chilean eastern boundary current region, often broaching the surface of 
the ocean in spite of predominantly equatonv'ard winds throughout the year. Figure 1, 
from Silva and Fonseca (1983), illustrates the relationship between the surface current 
and undercurrent at different times of the year. The PCU is overlain by an equatorward 
surface current south of 30°S, but broaches the surface and widens north of 30°S. The 
PCU is wider and more of a dominant feature in the north during the summer months, 
and becomes narrower during the winter. In the south, the PCU is dominated by the 
surface current, especially in the summer months, when winds are more meridional in 
orientation. Tchernia (1980) also acknowledges the existence of an undercurrent, de- 
scribing it as distinct and flowing beneath the countercurrent, with depths extending to 
several hundred meters. 

Until recently, the PCU had been relatively unexplored. Gunther (1936) dis- 
cussed the undercurrent that would eventually carry his name, but very little additional 
research was conducted on the PCU until 1961, when Wooster and Gilmartin (1961) 
measured it using parachute drogues. Their article sparked a new interest in the region, 
and soon several expeditions set out to study the waters adjacent to the South American 
coast. Between 1960 and 1982 the area was explored by twelve different expeditions 
(Fonseca, 1989). Results of the experiments revealed some consistent traits of the PCU; 

1. The PCU is present all year. 

2. It is easily distinguishable from the surface to about 600 m depth. 

3. Its core is at approximately 150-200 m depth. 

4. It is principally located between the coast and 200 km offshore. 

5. Average speeds are on the order of 25 cm/s. 

6. Transport is estimated to be 10 x 10* mV^. 

In a recent workshop on poleward flows along eastern ocean boundaries, 
Clarke (1989) concluded that, although many mechanisms may cause mean flows along 
eastern boundarv’ regions, only two have thus far been shown to have a significant im- 
pact on the creation of a poleward flowing undercurrent. The first of these is the gen- 
eration of poleward flow' by small-scale irregularities on the ocean bottom, in w'hich 
coastally trapped lee waves are formed by the irregularities under the influence of a mean 
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cquatorward flow. Poleward flows of this type were estimated by Clarke (1989) to be 
of such a small magnitude (~ 2.5 cm/s) that we will neglect them in this study. The other 
mechanism mentioned by Clark (1989) that can generate poleward flows of significant 
magnitudes (~ 25 cm/s) is wind stress, the role of which will be studied in this paper. 

B. CLIMATOLOGICAL WINDS 

The Chile Current System is influenced predominantly by equatorw^ard winds 
throughout the year. The direction, and to a large degree the magnitude of the winds 
are the result of a semi-permanent subtropical high pressure system that is similar in 
nature and behavior to its counterpart in the Northern Hemisphere, the North Pacific 
Subtropical High, as described in Nelson (1977). 

The South Pacific Subtropical High migrates meridionally with the seasons, reaching 
its northernmost extent during the austral winter, with its center located roughly at 
20® S. As the seasons progress, the high pressure system moves southward until its cen- 
ter reaches about 30®S during the summer months. 

Because of this migration, maximum values of wind stress and wind stress curl vary 

temporally at given locations. In Peru, for example, the maximum wind stress and wind 

stress curl values are seen during the austral winter months (Bakun and Nelson, 1991), 

# 

presumably because the South Pacific Subtropical High reaches its northernmost mi- 
gration at that time, and the resulting winds take on more of an alongshore character- 
istic. The middle of the Chilean coastline, in contrast, does not experience maximum 
wind stress and wind stress curl values until late spring and summer, corresponding to 
the upwelling season in that region (Fonseca, 1989). Maximum wind stress curl values 
for the entire region are seen during the austral summer months, near 33°S. 

Winds stress curl values are usually at a relative minimum in the region off the 
northwest coast of Chile, compared to the rest of the Chilean coastal waters. Bakun and 
Nelson (1991) argue that this curl minima (i.e., anticyclonic curl) might be caused by the 
coastal winds blowing into the coastal bend near Arica, rather than remaining parallel 
to the coast. As the winds blow onto the land, alongshore speeds would in effect go to 
zero, eliminating any wind stress curl-induced upwelling. Wooster (1970) discusses a 
cyclonic curl minimum in the same area, although he concentrates on the curl of the 
surface current, rather than on that of the wind. His discussion does, however, tend to 
support the observation that cyclonic curl values are at a relative minimum near Arica. 
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Figure 1. Surface current and undercurrent, from Silva and Fonseca (1983): Sim- 
plified representation of the Peru-Chile Undercurrent (shaded) between 
the coast of Cliile and ^ 290 km offshore. 
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III. MODEL DESCRIPTION 



A. MODEL EQUATIONS 

The numerical model used in this study was originally developed by Haney (1985), 
and later adapted by Batteen et al. (1989) for application to eastern boundary current 
regions. The following description of the model equations is based primarily on similar 
descriptions found in those papers. The multi-level, primitive equation (PE) model uses 
the hydrostatic, rigid hd, /?-plane, and Boussinesq approximations. The governing 
equations are: 



Equations of Motion : 



du —1 . /. . yA rr u . 



Po dx 



dz^ 

2 






dv -1 ^ d^V . , r ■ 

= 'Tn ~ V V -I- + (5^V. 



Po dy 



«5<i(v) 



oz 



w = — 



-H 



f du , dv s , 

( ^T” + T~ 
ox dy 



[1] 

[ 2 ] 

[3] 



Hydrostatic Equation: 

Equation of State: 
First Law: 




J_ 

H 



(*Q 




pgdc)dz 






p = Po(i-«(r-7o)) 

02 



[4] 

[5] 

[ 6 ] 



Important features of the above equations include: 

• A right-handed cartesian coordinate system is used in the above equations, with x 
pointing in the easterly direction, y pointing alongshore in a northerly direction, 
and z pointing upward. Accordingly, u is velocity in the x direction, and is positive 
when pointing toward the east; v is positive to the north, and w is positive upward. 

• T is temperature. 

• p is density. 

• /is the Coriolis parameter. 

• p' is the departure from the vertically averaged pressure (which is assumed to be 
zero). Because the mean pressure is assumed to be zero, the barotropic mode is 
considered to be ignored. 

• £ is simply a constant of integration used in Equations (3) and (4). 
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• Equation (5) assumes p to be a function of temperature only. Any dependence on 
salinity is ignored, 

• The terms 6J^u), <5<Xv) and bJJ) represent the vertical turbulent mixing of zonal and 
meridional momentum and heat, respectively, due to surface layer processes, by a 
dynamic adjustment mechanism. This adjustment maintains the dynamic stability 
in the water column, in the sense of lighter water overlying denser water (Haney 
et ai, 1978). 

The boundary conditions at the top (z = 0) of the model ocean are : 





[7] 


dv T 

dz Po 


[8] 


Kh-^ = -Qs 


[9] 


IV = 0 . 


[10] 


In Equation (8), t is the surface wind stress. In Equation (9), 
of longwave radiation. 


is the net upward flux 


The bottom (2 = —H) boundary conditions are as follow: 


# 


-|^ = C^(w^ + vy\u cos y - V sin y) 


[11] 


Kfn cos y + M sin y) 


[12] 


PT 

Kh-^ = 0 
02 


[13] 


w = 0 . 


[14] 



In the above equations, Q is a bottom drag coefficient and y = 10® is the geostrophic 

inflow angle (Weatherly, 1972). Table 1 provides a list of other symbols found in the 

model equations, as well as values of constants used throughout this study. 

1 

B. MODEL DOMAIN AND RESOLUTION 

The domain of both experiments in this study consists of a rectangular area, with a 
meridional length of 1280 km, and a zonal width of 576 km. The eastern side of the box 
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Table 1. VALUES OF CONSTANTS USED IN THE MODEL 





Value 


Name 


a 


2 n da\r^ 


Earth rotation rate 


c 


0.958 cal gm~^ ®/f“* 


Specific heat of sea water 




1.22 X 10-' 


Bottom drag coefficient 


c 


1.3 X 10-' 


Atmospheric drag coefficient 


L 


595 cal gm-' 


Latent heat of sea water 


To 


21S.2°K 


Constant reference temperature 


Pc 


1.23 X 10-' gm cm~^ 


Density of air 


Po 


1,0276 gm enr^ 


Density of sea water at 


a 


2.01 X 10-^ ®A'-' 


Thermal expansion coefficient 


K 


10 


Number of levels in vertical 


AA' 


9.0 X lO' cm 


Zonal grid spacing 


^Y 


2.0 X 10* cm 


Meridional grid spacing 


D or H 


4.5 X 10' cm 


Total ocean depth 


00 


22.0® S 


Latitude of northern boundary 


4>c 


34.0® S 


Latitude of southern boundary 




-72® W 


Longitude of eastern boundary 




- 78® W 


Longitude of western boundary 


At 


800 s 


Time step 


g 


980 cm 5 -' 


Acceleration of gravity 


Am 


2.0 X lO’’ c/7p 5 -' 


Biharmonic momentum dilfusion coefi'. 


A„ 


2.0 X lO’’ enf 5 -' 


Biharmonic heat diffusion coefficient 


Km 


0.5 enP 5 -* 


Vertical eddy viscosity 


K, 


0.5 cm^ 5 -' 


Vertical eddy conductivity 



is assumed to touch the South American coast (Figure 2). In addition, for modeling 
purposes, the coastline is assumed to run exactly north-south; consequently, the domain 
is similarly oriented. The actual coastline is relatively straight, but runs in a slight 
southvvcst to northeast direction. At 34® S the western edge of the coast is approxi- 
mately 72® W; at 22® S, the coast is located at about 70® W. 
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C. FINITE DIFFERENCE SCHEME 



The numerical model used in this study is a staggered grid model, using the Arakawa 
and Lamb (1977) B-scheme. Batteen and Han (1981) have shown that this scheme is 
appropriate when grid spacing is approximately on the same order as, or less than, the 
Rossby radius of deformation, which meets the criteria of this study. 

In the vertical, the model has ten layers, with constant z-levels of 13, 46, 98, 182, 
316, 529, 870, 1416, 2283, and 3656 m. Consistent with Haney (1974), this spacing 
scheme concentrates more on the upper, dynamically active part of the ocean, above the 
thermocline. 

D. HEAT AND MOMENTUM DIFFUSION 

Because this study is designed to investigate the generation of mesoscale eddies 
along the coast of Chile, the type of diffusion used in the model equations becomes an 
important issue. Laplacian lateral heat diffusion tends to suppress the formation of 
mesoscale eddies by diminishing baroclinic signals at that scale, thus making it more 
difficult for baroclinic instabilities to overcome diffusive damping (Holland and Batteen, 
1986). Biharmonic diffusion, on the other hand, tends to allow the formation of 
mesoscale instabilities, because it tends to damp at scales smaller than mesoscale eddies 
(Batteen et aL, 1989). For these reasons, biharmonic lateral heat and momentum dif- 
fusion are used in this paper. Coefficients are listed in Table 1. 

E. SURFACE THERMAL FORCING 

To isolate the effects of wind forcing in these process-oriented studies, the net heat 
flux at the sea surface should initially be zero. In this way, any heat flux that is gener- 
ated will have been generated by an increase or decrease of sea surface temperature re- 
sulting from wind forcing effects (Batteen et ai, 1989). To accomplish this, surface 
generated heat flux {w'T') is set equal to (5^) , which is defined as the solar radiation at 
the sea surface (Haney et ai, 1978). Simply stated: 
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So-(wT') = 0. 



[15] 



These terms and their components are in turn broken down into the following math- 
ematical definitions: 



^0 — ^0 I PqC C^6] 

{^)o = {Qb + Qs + Qe)IPoC [17] 

So = 0.95(0.74 - 0.6/;)S^ [18] 

S^ = 3.642 X 10”^ [19] 

<?e = 0.985(7(7*5)'* X [0.39 - 5.0(c^)*^^](l - 0.6«^) [20] 

Qs = PACHCp\y\{Ts - T,) [21] 

Q, = I V I (0.622/p J(C5 - cj [22] 



The variables used in these equations are defined in Table 2. Values for Ts,n,p^,e^, 
specific and relative humidity (used to compute ) are all taken from Esbensen and 
Kushnir (1981), while solar insolation values are taken from List (1984). To estimate 
annual means, the values were averaged over twelve months. Using these values, the 

ambient air temperature that results in zero heat flux at the surface is 290.13® K. 

F. BOUNDARY CONDITIONS 

The rectangular box surrounding the model domain is bordered on the northern, 
southern and western edges by open ocean. The eastern edge is assumed to border the 
western coast of Chile; it is a no-slip, straight vertical wall extending from the surface 
to 4500 m depth. At this eastern boundary, both tangential (y direction) and normal (x 
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Table 2. VARIABLES USED IN COMPUTING NET HEAT FLUX 





Value 


Description 


So 


Equation 16 


Solar radiation at sea surface 


w'T' 


Equation 17 


Surface generated heat flux 


So 


Equation 18 


Downward flux of solar rad 


s: 


Equation 19 


Flux of solar radiation at top of 
atmosphere 


Qb 


Equation 20 


Net upward longwave rad 


Qs 


Equation 21 


Sensible heat flux 


Qb 


Equation 22 


Latent heat flux 


T 


1 dynejan^ 


Wind stress 


C 


0.958 cal r* 


Spec heat of sea water 


Q 


0.24 cal g-* 


Spec heat of dry air 


n 


0.63 


Fractional cloud cover 




5% 


Reflectivity of sea surface 




26 % 


Reduc of insol by atmos absorp 




1.757 X 10-^ 


Heat exchange coeflicient 


|V| 


830 cm/s 


Wind speed 


T, 


16.72° C 


Sea surface temp 


T. 


290.13° K 


Ambient air temp 


Pa 


.00123 g cm~^ 


Air density at sea level 


a 


1.35 X 10-‘2 


Stefan-Boltzmann constant 




20.8 mb 


Saturation vapor press at 


e. 


16.49 mb 


Atmos vapor pressure 


Pa 


1013.25 m^ 


Atmos press at sea level 



direction) components of velocity go to zero. A flat bottom at 4500 m is assumed 
throughout the entire domain. Both of the above assumptions can be justified for two 
reasons. First, because this is a process-oriented study, the inclusion of a rough 
coastline or bottom would induce additional shears and curl values (compared to those 
induced by a flat bottom and straight vertical coastline as in Batteen et al. (1989) and 
Fonseca (1989)), with the result that the particular process being investigated (i.e., wind 
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forcing) would not be isolated. Second, the Chilean coastline is relatively straight and 
smooth, the continental shelf is relatively narrow, and the continental slope is steep, 
dropping from sea level to over 1000 fathoms (~ 2000 m) in a distance of less than 
25 km from the coast (Ziegler et ai, 1957). Therefore, the use of a straight, vertical wall 
is not too unreasonable an approximation for a process-oriented study designed to iso- 
late the effects of wind forcing on the Chilean Current System. 

G. INITIAL CONDITIONS 

Climatological vertical and horizontal temperature profiles, derived from Levitus 
(1982), show that, over the study domain, the average seasonal surface temperatures 
vary between a summertime maximum of roughly 22® C at ^ 22® S, to a wintertime 
minimum of 13® C at ^ 34® S. The average annual sea surface temperature for the area 
is 1 7® C. Although a relatively weak north-south temperature gradient is evident at all 
depths, seasonal variations are negligible below 100 m depth. 

A temperature profile based on Levitus (1982) was used in both experiments to 
represent the initial temperature conditions for the model domain. The method for in- 
corporating the values of the profile into the model was similar to that used in Batteen 
and Rutherford (1990). Initially, mean annual temperatures were taken from Levitus 
(1982) for applicable 5° squares in or adjacent to the model domain, then interpolated 
in an east-west direction to match the middle of the domain. The temperature profile 
at each point was then interpolated vertically to correspond to the midpoints of the ten 
layers used in the model. These computations resulted in two "smoothed" temperature 
profiles, one representing the northern part of the domain (centered at 17.5® S), and the 
other representing the southern region (centered at 37.5® S ). These two profiles were 
interpolated once again, this time in a north-south direction, to obtain an annual tem- 
perature profile, centered at 28° S. 
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As in Batteen et al. (1989, 1991) an exponential curve was then used to approximate 
the climatological temperature profile that had been computed from Levitus (1982). The 
form of the equation used was: 

T{z) = Ts + ATe"\ [23] 

The approximation uses a length scale of 450 m, and assumes Ts — 2°C to be the tem- 
perature at great depth. AT = 15° C is the change of temperature between the bottom 
and top of the profile. Figure 3 shows that the resulting temperature profile fits the 
mean annual interpolated profile. This exponential profile was used to provide initial 
values at each of the ten levels in the model. 

H. EXPERIMENTAL CONDITIONS 

1. Experiment One; Constant Wind Stress 

In the first experiment, a constant equatonv’ard wind stress of 1 dynelcm^ 
( 830 cm/s) is applied to all gridpoints within the domain, except for a band at the 
northern and southern ends, where no wind forcing is applied. McCreary (1981) dem- 
onstrated that this technique allowed for the propagation of coastal Kelvin waves, which 
in turn allowed for the creation of an undercurrent and a coastal jet (Batteen et al., 
1989). Because the winds are predominantly equatoru^ard all year off Chile, the model 
is run for 240 days, which is a sufficient time for currents and eddies to develop in the 
model. 

2. Experiment Two: Variable Climatological Winds 

In Experiment 2 , daily climatological winds are used to simulate the temporal 
and spatial wind variations that have been observed off the coast of Chile. Data, in the 
form of climatological wind stress fields, are taken from Bakun and Nelson (1991). The 
wind stress fields were derived from all available ship wind reports that have been col- 

I 

lected in the region between 1850-1979, after eliminating physically unrealistiV- reports. 
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Figures 4 and 5, taken from Bakun and Nelson (1991), show the wind stress and wind 
stress curl fields, respectively, that are used. The wind stress fields are grouped into 
one-degree square regions (their midpoints being on the intersections between integer 
latitudes and longitudes), and a monthly mean value is computed for each region. 

Because wind speed values (vice wind stress values) are used in the model, wind 
speed values were obtained using the following equation: 

""VS- 

where t is wind stress, is air density, is the atmospheric drag coefficent, and U is 
wind speed. The resulting wind values were then smoothed using a nonlinear 3X3 
median smoothing technique described in Rabiner et al. (1975), and used successfully in 
Batteen et al. (1991). This technique was shown in Rabiner et al. (1975) to preserve the 
significant features in a field that can otherwise be hidden or "smoothed-out" when a 
linear filter is used. After smoothing, the resulting wind values for the 
n” latitude by 6® longitude region being studied were interpolated onto a 65 X 65 grid, 
basically covering the same geographical area. The easternmost gridpoints were con- 

t 

sidered to be touching the continent. The resulting twelve monthly w ind fields were then 
linearly interpolated into 365 daily fields, which were used to simulate one year's worth 
of daily, climatological wind values. In their final form, the data consisted of 365 grids 
(65 X 65), each gridpoint having a unique interpolated value. In the model, the appro- 
priate day's values are read at the beginning of each integer day. An example of the 
wind fields can be seen in Figure 6, which show's the annual migration of maximum 
winds from the south in the austral summer, to the north in the austral winter months, 
the effects of which are discussed in Chapter V. 

It should be kept in mind that the daily wind values are a linear interpolation 
of monthly climatological winds, and not daily climatological values. As a result, the 
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values should be considered as linear extensions of month-to-month trends, divided into 
daily values to allow for a smooth transition from one monthly mean to the next. 
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Figure 2. Model Domain: Area of study is a 516 km by 1280 A:m box off the 
coast of Chile. 



16 





Figure 3. Temperature proFiie used in both experiments: A smooth exponential 
curv'e (solid line) was used to approximate the climatological temper- 
ature profile (solid line with small circles) computed from 
Levitus (1982). 
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Figure 4. Wind stress fields from Bakun and Nelson (1991): Vectors are sealed in 
units of dynelcm^ . Ellipses show areas of large standard errors (dis- 
cussed in Bakun and Nelson, 1991). 
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Figure 5. Wind stress curl fields from Bakun and Nelson (1991): Units and con- 

tour interval arc lO'* dynelcm ^ . Regions of anticyclonic curl are shaded. 
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Figure 6. Wind speed vectors with magnitude contours: Vectors are scaled in units 

of cm/s. Contours are in hundreds of cm/s. Day 1 corresponds to 1 Jan- 
uary. 
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IV. ANALYSIS TECHNIQUES 



The following energy and spectral techniques, adapted from Batteen and Rutherford 
(1990), will be used in the following chapter to analyze the generation and stability of 
the Peru-Chile Current and eddies. 

A. ENERGY ANALYSIS 

An energy analysis based on that of Han (1975) and Semtner and Mintz (1977) is 
made to gain a better understanding of the energy transfers in the unstable flow. The 
energy calculations are presented using the Semtner and Mintz (1977) notation: 



and presented in a time scries plot. After reaching a quasi-steady state in which the total 
kinetic energy is nearly constant, the time mean and time eddy kinetic energy are calcu- 
lated by: 



( ) time average 

( )' time deviation 

( ) horizontal space average 

( )* horizontal space deviation. 



The kinetic energy (K) is calculated by: 




[25] 




[26] 




[ 27 ] 
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Available potential energy (P) is calculated by: 



P = ag 



T'Mirr 



[ 28 ] 



The temporal mean and eddy available potential energy are calculated by: 



P = or g 






[29] 









[30] 



The transfers between the energy types are defined, after Semtner and Mintz (1977), 



by: 



{/f ^ />) = _ «g[ru;] 


[31] 


{/>' K') = 


[32] 


[K K'} = v.^V . v'v' + 


[33] 


{P^P') = ag 


[34] 



The model output consists of velocity components and temperature at daily intervals 
at each gridpoint and neither the vertical velocity nor any advection terms calculated by 
the model are stored. Because the calculation of the energy transfers requires both ver- 
tical velocity and numerous advection terms, these were recalculated in the same manner 
as in the model but using the stored values of u, v and T as the input data. Although 
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the energy transfers are recalculated from model output, they are consistent with the 
initial calculations of vertical velocity and advection terms in the model. 

Semtner and Mintz (1977) applied their energy transfer analysis to currents which 
had become unstable, generated eddies and then reached a quasi-steady state. In this 
study, the quasi-steady energetic state prior to eddy generation is examined, and the 
energy transfers are used to argue for the instability mechanism which leads to the initial 
eddy generation in each case. 

B. SPECTRAL ANALYSIS 

A spectral analysis technique is used to estimate the dominant wavelengths at which 
eddy growih occurs. The alongshore wavenumber is selected based on an assumption 
of a meridional anisotropic preference for the eddy development. The model configura- 
tion of 65 alongshore temperature gridpoints and 64 velocity points is ideal in allowing 
data over the entire domain to be used to compute the discrete Fourier transform based 
on 64 points. Leakage due to the finite length of the series is reduced using a Hanning 
window. The loss of variance due to the w'indowing is then compensated by scaling the 
calculated one-sided spectral density appropriately. 



r 
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V. RESULTS OF THE WIND FORCING EXPERIMENTS 



Experiment 1 studies the effects of steady, upvvelling-favorable winds in the Chile 
eastern boundary current region. Experiment 2 investigates the effects of seasonally 
varying winds in the same region, using climatological monthly values that have been 
interpolated to daily values. 

A. EXPERIMENT 1 

In Experiment 1, a constant uind stress was applied to the interior of the domain 
for 240 days. As in Batteen et al. (1991), inertial oscillations of near-surface ocean cur- 
rents developed initially. After a few days of model simulation time, the oscillations were 
damped, leaving quasi-steady offshore Ekman transport, as shown by the surface veloc- 
ity vectors in Figure 7A. Figures 7B, 7C and 7D indicate that an equatoru'ard-flowing, 
upwelling-favorable surface current had been generated. A poleward-flowing undercur- 
rent (not shown) had also developed. The surface current flowed at ~ 20 cm/5, and ex- 
tended ~ 50 km offshore, while the undercurrent, centered along the coast at ~ 300 m 
depth, had weaker velocities of ~ 4 cm/5. 

The first evidence of developing eddies occurred by day 50, as evidenced by the ap- 
pearance of perturbations in the zonal velocity fields (Figure 8). (Perturbations are 
more clearly seen in these types of presentations than in meridional velocity plots, be- 
cause the latter tend to be dominated by the currents, which are predominantly in the 
meridional direction.) The time sequence of surface zonal velocity shows the develop- 
ment and growth of eddies from days 50-160. The baroclinic and barotropic transfers 
(based on Equations (32) and (33), respectively), are given in Figure 9, and show that 
both types of instabilities contribute to the development of eddies in the region. 
Barotropic instability is evident near the offshore region of the coastal jet and near the 
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coastline because of the horizontal shear in the region. Baroclinic instability, which 
draws its energy from the vertical shear in the mean flow (Haltiner and Williams, 1980; 
Batteen et ai, 1991), was evident in the interior of the coastal region. 

By day 60, meanders were evident at y ~ 250 to 500 km near the coast in the velocity, 
dynamic height and temperature fields, as seen in Figures lOA-lOD. As in the )?-plane 
experiment of Batteen et ai (1989), eddies first formed in the poleward end of the region, 
as expected. Figures 10 and 11 show that the dominant wavelength at day 60 is 
~ 100 km. 

Figure 12 shows three meridional velocity cross-sections averaged over days 50 - 60. 
The three cross-sections correspond to areas near the equatorward, middle, and 
poleward parts of the model domain, respectively. Figure 12 shows that the surface 
current has a width of ~ 100 km, and an average velocity of ~ 18 cm/s, while the under- 
current has a core velocity of 10 cm/s at ~250m depth. The undercurrent weakens 
towards the equator and intensifies near its poleward end. 

Figure 13 shows the fields of surface velocity vectors, temperature, velocity con- 
tours, and dynamic heights (relative to 2000 m) at day 120. The perturbations and 
meanders seen earlier in the southern part of the domain (see Figure 7) have developed 
into anticyclonic and cyclonic eddies. For example, a cyclonic eddy is discernible at 
y ~ 256 km with a corresponding cold core temperature of ~ 14.5° C. Figure 14 shows 
that, although both barotropic and baroclinic instability processes contribute to the 
formation and development of eddies, barotropic energy transfer is dominant, due to the 
larger energy transfer. Figure 15 indicates that the surface current has intensified to 
> 20 cmjs near the poleward end of the domain. The undercurrent maintains a core 
velocity of ~ 5 — 10 cmjs , and is better developed at ~ 620 km than towards the pole. 

Closed eddies were fully formed by day 160 (Figures 16A and 16D). While average 
values of the coastal surface current were ~ 30 cmjs (not shown), a maximum value of 
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~ 40 cm/s was reached by day 160 (Figure 16C). A cold, offshore-flowing filament ex- 
tends to ~ 350 /cm offshore (Figure 16B). Figure 17 shows that the current extends to 
~ 150 — 200 m depth in the equatorward end of the domain, while the undercurrent core 
is at ~ 350 m and centered at 10 km offshore (except at the poleward end, where the 
core is at the coast). Figure 18 shows that the dominant wavelength w^as ~ 250 km. 

The maximum velocity of the undercurrent was reached by day 210, when it reached 
a magnitude of ~ 25 cm/s (Figure 19). The core moves from 10 km (Figure 19A) to 
~ 40 Am (Figure 19C) offshore from north to south. What appeared to be a surface 
countercurrent in Figure 19 is most likely a manifestation of a cyclonic eddy that had 
formed in that area, with its poleward-flowing eastern side dominating the 
equatorward-flowing mean surface flow. This is supported by Figure 20D, which shows 
an upwelling-favorable low height center in the same region {y 300 Am), accompanied 
by a corresponding cold core center. Figures 20A and 20C both show poleward com- 
ponents in the vicinity of the eastern edge of the cyclonic eddy. Also evident are cold, 
oflfshore-flowing filaments extending as far as ~ 400 Am offshore, as seen in Figure 20B. 

B. EXPERIMENT 2 

In Experiment 2, daily-var\ing wind forcing, corresponding to interpolated monthly 
climatological values taken from Bakun and Nelson (1991), were used in the model. The 
experiment ran from day 1 (corresponding to January 1) to day 365. In the early part 
of the experiment, which corresponded to the middle of the austral summer, winds were 
generally meridional throughout the domain (see Figure 6), and the South Pacific 
anticyclonic gyre was at its southernmost extent. Figure 5, taken from Bakun and 
Nelson (1991), also shows a relative wind stress curl minima in the northeast corner of 
the domain during the summer months. This would imply that upwelling would be en- 
hanced in the southern portion of the domain, w^here cyclonic wind stress curl values 
were higher. As was the case in Experiment 1, an equatorw'ard surface current was es- 
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tablished by day 20, extending 80 km ofT the coast and reaching a magnitude of 
> 20 cm/5. An undercurrent first appeared at day 20 in the southern extreme of the 
region, but by day 40 was present all along the coast, with a core 10 km offshore, at 
350 m depth (see Figure 21), 

A series of temperature cross-sections at day 50 (Figure 22) shows upward bending 
of isotherms near the coast, which is consistent with a more intense upwelling regime 
in the southern than in the northern part of the region. A planar view of the surface 
temperature at day 50 (Figure 23B) also shows a wider sw^ath of cooler, upwelled water 
in the south, and a narrower band in the north. Surface velocity vector, zonal velocity, 
and dynamic heights in Figures 23A, 23C and 23D, respectively, show that, in contrast 
to Experiment 1, eddies had not yet formed by day 50, and velocities tended to be 
smaller, perhaps because of the weaker winds (generally less than 800 cmis) throughout 
most of the region. 

The first evidence of perturbations w^as seen by day 80 (Figure 24), Meanders occur 
~ 100 Am to 200 km further to the north than their counterparts in Experiment 1. The 
reason for the difference between the tw^o experiments was due to the different wind 
forcings used (i.e., uniform versus seasonally varying winds), because all other conditions 
remained the same. Spectral density analysis reveals a smaller w'avenumber than had 

f 

been present in Experiment 1, corresponding to a w^avelength of ~ 85 Am (Figure 25), 
(also seen in Figure 24, when eddies are being generated). Figure 26 shows that by day 
80 a w'eak, poleward flowing undercurrent is only evident in the middle and southern 
portions of the region, centered w ithin 20 Am of the coast, at a depth of 300 m. Also 
shown is a surface current extending ~ 80 Am off the coast, to a depth of ~ 150 m . 

By day 140, the South Pacific High had migrated to the north, allowing a westerly, 
diffluent wind pattern to develop in the southern part of the domain (Figure 6), while 
winds in the north remained meridionally oriented. A time sequence of dynarhdc height 
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fields (Figure 27) shows that eddies were well developed by day 140, with larger eddies 
forming in the north, and smaller eddies in the south. As discussed in Batteen and 
Rutherford (1990), in a ^-plane regime, barotropic instability tends to support eddies 
with radii greater than the Rossby radius of deformation, whereas baroclinic instability 
tends to support growth on the same order as the Rossby radius of deformation. As 
Figure 28 indicates, barotropic instability is prevalent through most of the coastal re- 
gions, w'hile baroclinic instability tended to dominate only in the extreme south, below 
the region of diffluent winds. Because this effect was not seen in Experiment 1, the 
varving winds must have been the cause of the different sized eddies. By day 140 the 
swath of cooler, upwelled temperatures was no longer wider in the south (Figure 29). 
Instead, the widest band of cool temperatures was found in the north, once again likely 
due to the diffluent winds in the south, which tended to blow towards the coast rather 
than along it (Figure 6). In Figure 29, the time sequence of sea surface temperature 
contours also show's the development of cold filaments from days 80 to 140. 

By day 250, the winds had transitioned back into primarily a meridional regime 
(Figure 6). The time-averaged surface current had increased to > 50 cmjs along the 

t 

middle of the coast by day 250 (Figure 30), where wind magnitudes were greatest (Fig- 
ure 6). The dominant wavelength (Figure 31) had increased to ~ 250 Am, which cor- 
responded to the large eddies in the north. Smaller eddies were still evident in the south, 
and, as indicated in Figure 31, had become nearly as energetically significant as the 
eddies at longer wavelengths. The smaller scale eddy generation occurred during a pe- 
riod of primarily meridional, equatorw^ard winds. 

Figure 32 show's a time sequence of surface dynamic height fields from days 250 to 
280, w'hen a cyclonic, cold-core eddy at ~ 700 km breaks off and transits to the west. 
Figure 33 shows a time sequence of dynamic height fields towards the end of the year, 
W'hen an anticyclonic, warm core eddy at ~ 768 km travels eastward, into the region 
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of lower dynamic heights, possibly in response to a more meridionally oriented wind 
field. 

The equaton^’ard surface current, which reached a maximum speed of ~ 50 cmis by 
late austral spring (when winds had returned to a meridional orientation throughout the 
region), remained relatively constant throughout the duration of the experiment. The 
average magnitude of the undercurrent (averaged over 30 days to filter out eddies) was 
~ 15 cm/s, at the end of the model simulation, when the South Pacific High had shifted 
back to the southern part of the domain. The largest instantaneous velocity of the 
undercurrent (not shown) was ~ 32 cm/s , occurring in the southern part of the region 
at day 360. 
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Figure 7. Surface current vectors, temperature, meridional velocity and dynamic 
height at day 50 of Experiment 1: Velocities are in cm/s, T is in ®C, 
heights are in cm relative to 2000 m. Dashed lines are negative height 
values. 
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Figure 8. U (surface zonal velocity) for Experiment 
Dashed lines denote westward velocities. 



Velocities are in cm/s. 
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BAROCLINIC TRANSFER 
^^sec^'1 
Chile steady wind 
Days 40—50: Layers 1 — 5 




BAROTROPIC TRANSFER 
(erqs.cm'^.sec'^') 

Chile Wind830 
Doys 40 —50: Layers 1 - 5 





Figure 9. Baroclinic and barotropic instability at day 50 of Experiment 1. 



32 




C at doy 60 
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T ot day 60 depth 0 




P at day 60 depth 0 




Figure 10. Surface current vectors, temperature, meridional velocity and dynamic 
height at day 60 of Experiment 1: Velocities are in cni/s, T is in ®C. 

heights are in cm relative to 2000 m. Dashed lines denote negative 
values. 
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Day 6Ck 




Figure 1 1. Spectral density at day 60 of Experiment 1: Dominant wavelength is 

~ 100 km, where X = wavenumber"* . 
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Figure 12. Zonal cross-sections of north-south velocity (V) at day 60 of Exper- 
iment 1: Velocities are in cm‘,s. Southward flowing currents are 

dashed. 



35 




A 



C 



C at day 120 





V at day 120 depth 0 





Figure 13. Surface current vectors, temperature, meridional velocity and dynamic 
height at day 120 of Experiment 1: Velocities are in cm/s, T is in ®C, 
heights are in cm relative to 2000 m. Dashed lines denote negative 
values. 
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Figure 14, 



Baroclinic and barotropic instability at day 120 of Experiment 1 





Figure 15. Zonal cross-sections of north-south velocity (V) at day 120 of Exper- 
iment 1; Velocities are in cm',s. Southward flowing currents are 
dashed. 
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Figure 16. Surface current vectors, temperature, meridional velocity and dynamic 
height at day 160 of Experiment 1: Velocities are in cm/s, T is in "C, 
heights are in cm relative to 2000 m. Dashed lines are southward ve- 
locities (C) or negative values (D). 
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Figure 17. Zonal cross-sections of north-south velocity (V) at day 160 of Exper- 
iment 1: Velocities are in crnis. Southward flowing currents are 

dashed. 
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SPECTRAL DENSITY 
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Figure 18. 



Spectral density at day 160 of Experiment 1: 



Dominant wavelength 



is ~ 250 km. 
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Figure 19. Zonal cross-sections of north-south velocity (V) at day 160 of Exper- 
iment 1: Velocities are in cm/s. Southward flowing currents are 
dashed. 
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Figure 20. Surface current vectors, temperature, meridional velocity and dynamic 
height at day 210 of Experiment 1; Velocities are in cm/s, T is in ®C, 
heights are in cm relative to 2000 m. Dashed lines are southward ve- 
locities (C) or negative values (D). 
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Figure 21. Zonal cross-sections of north-south velocit)’ at day 40 of Exper- 
iment 2: Velocities are in cm/s. Southward flowing currents are 

dashed. 
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Figure 22. Zonal cross-sections of north-south temperature (T) at day 50: Tem- 
peratures are in ®C. 
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Figure 23. Surface current vectors, temperature, zonal velocity and dynamic height 
at day 50 of Experiment 2; Velocities are in cm/s, T is in ®C, heights 
are in cm relative to 2000 m. Dashed lines are westward velocities (C) 
or negative values (D). 
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Figure 24. 



Surface zonal velocity (U) for Experiment 2: 



Velocities are in cm/s. 



Dashed lines are westward velocities. 
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SPECTRAL DENSrTY 

Experiment I (STEADY Wind) 



Day 80: 




Figure 25. Spectral Density at day 80 of Experiment 2: Dominant wavelength 

is ^ 250 km, with a secondary wavelength evident at ~ 85 km. 
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Figure 26. Zonal cross-sections of north-south velocity (V) at day 80 of Exper- 
iment 2: Velocities are in cmjs. Southward flowing currents are 
dashed. 



49 




A 



P at day 80 depth 0 




Distance (km) 




contour interval = 2.0 



contour interval = 2.0 




Distance (km) 



D 



P at doy 140 depth 0 

1280 



1024 



768 J 
O 

o 

c 

o 

-512 

O 



256 



0 

576 384 192 0 

Distance (km) 




contour interval = 2.0 



contour intervol = 2.0 



Figure 27. Dynamic height, days 80-140 of Experiment 2: Heights are in cm rel- 

ative to 2000 m. Dashed lines denote negative height values. 
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Figure 28. Baroclinic and barotropic instability at day 140 of Experiment 2 . 
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Figure 29. Surface temperature contours for Experiment 2: Temperatures are 
in ®C. Temperature bands in the south are narrower during the austral 
fall. 
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Figure 30. Zonal cross-sections of north- south velocity (V) at day 250 of Exper- 
iment 2: Velocities are in cni/s. Southward flowing currents are 

dashed. 
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SPECTRAL DENSrTY 

Experiment I (STEADY Wind) 




Figure 31. Spectral density at day 250 of Experiment 2: Dominant wavelength 

is ~ 250 km, with a secondary w'avelcngth evident at 85 km. 
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Figure 33. Dynamic height, days 290-350 of Experiment 2: Heights are in cm 

relative to 2000 m. Dashed lines denote negative height values. 
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Figure 32. Dynamic height, days 250-280 of Experiment 2: Heights are in cm 

relative to 2000 m. Dashed lines denote negative height values. 
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VI. SUMMARY AND RECOMMENDATIONS 



A. SUMMARY 

This study used a high-resolution, primitive equation numerical model to study the 
effects of steady and daily vaiy'ing winds on the Chile Current System. The model do- 
main was a rectangular region extending 6° off the coast, between 22” and 34” S . 

In Experiment 1, a constant wind stress was applied to the model domain, excluding 
the northern and southern open boundaries, for a period of 240 days. The objective was 
to simulate the generally equatorvvard-flowing winds that are seen in that part of the 
world, for most of the year. Both an equatorward-flowing surface current and a 
poleward- flowing undercurrent developed. Average surface equatorward and subsurface 
poleward velocities were ~ 25 — 30 anjs. Eddies formed in the southern part of the do- 
main first, then spread north and south, eventually dominating the coastal region. The 
dominant wavelength for the eddies was ~ 250 km. 

In Experiment 2, temporally and spatially interpolated values of monthly 
climatological winds were used as the forcing for the model. Unlike Experiment 1, the 
varying winds allowed the following processes to be studied. First, it allowed the study 
of the effects of seasonal wind variations on the behavior of the current system. Second, 
it allowed the influence of wind stress curl on coastal dynamics to be investigated. Third, 
it introduced the zonal component of the wind stress, allowing the study of the oceanic 
response to this additional component of stress. 

Initially, when winds were generally equatorward, the two models were similar, with 
the development of coastal currents and upwelling along the coast. As the seasonal re- 
gime of the winds changed, differences, as expected, were discernible. For example, the 
poleward-flowing undercurrent in Experiment 2 was slower to develop, especially in the 
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northern parts of the region. The surface equatorward current reached a maximum 
speed of > 50 cm/s in Experiment 2 compared to a maximum speed of ~ 30 cm/s in Ex- 
periment 1. The average undercurrent in Experiment 2 was ~ IScmls, which is slower 
than in Experiment 1 (which had velocities of ~ 25 cm/s ). In Experiment 2, the under- 
current was centered at 300 m depth, compared to ~ 250 m in Experiment 1. In Ex- 
periment 2, instabilities did not appear until day 80, a full 30 days later than in 
Experiment 1. Baroclinic instability played a smaller role (than in Experiment 1) in the 
creation of eddies, except in the extreme southern part of the region. Once eddies did 
form, however, their sizes varied with latitude and time, tending to grow larger in the 
north where the winds were usually mostly meridional, and virtually disappearing in the 
south, when the alongshore wind stress was at a minimum. These southern eddies began 
to grow again, however, once stronger winds returned to the region. In the north, the 
dominant wavelength was ~ 250 km, while in the south it was ~ 85 Am. 

B. RECOMMENDATIONS 

Of critical importance to any process-oriented study is the comparison of model 
output with available observations. Without such comparisons, model studies are difii- 
cult to verify, and become speculative in nature. As mentioned by several authors (e.g., 
Fonseca, 1989; Bakun and Nelson, 1991), the waters olT the coast of Chile are 
notoriously lacking in both wind and current data, and until more observational studies 
are conducted, the Chile Current System will remain somewhat of a mysteiy. Johnson 
(1989) suggests that satellite imagery of the South American coastal waters may provide 
the most readily obtainable basis for increasing the number of available observations. 
Much of the model output from numerical simulations could then be compared directly 
to the imagery. Figure 34A shows a satellite image from Fonseca (1989) of a cold fila- 
ment extending out '>-40 Am off the coast of Chile, near 33° S, while Figure 34B shows 
the surface temperature field at day 290 of the Experiment 2 model run, in which cold 
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filaments can be seen extending to the west. The 14° contour at several locations 
protrudes 40-60 km from the coast, a distance very similar to that seen in satellite im- 
agery. Comparison of the two figures suggests that the model results may have some 
validity, although many more comparisons would need to be made before such a claim 
could be made with any certainty. This merging of model output and real-time satellite 
verification may well be the next logical step in understanding the Chilean coastal wa- 
ters, and is a recommended area of future study, as such observations (e. g., GEOSAT, 
SEASAT) become available. 

The knowledge gained in this investigation can benefit Naval Oceanography by 
providing a basis for prediction not only in the Chilean Current System, but also in other 
dynamically similar Eastern Boundary Current regions. An understanding of the for- 
mation and evolution of eddy and jet structures is an essential step in the development 
of a real-time analysis and prediction capability for these important features (Robinson 
et ai, 1984). 
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Figure 34. Temperature contours from (A) satellite and (B) model 
output: (A) Satellite-derived temperatures from Fonseca (1989); 
(B) Temperature contours from day 290 of Experiment 2. All tem- 
peratures are in ®C. The inside box at y ^ 256 km in (B) corresponds 
to the size of the region shown in (A). 
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